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Abstract: QCD analysis of ^^(a;,^^) is revisited. It is emphasized that the presence of the 
inhomogeneous term in the evolution equations for quark distribution functions of the photon 
implies important difference in the way factorization mechanism works in photon-hadron and 
photon-photon collisions as compared to the hadronic ones. Moreover, a careful definitions 
of the very concepts of the "leading order" and "next-to-leading order" QCD analysis of F2 
are needed in order to separate genuine QCD effects from those of pure QED origin. After 
presenting such definitions, I show that all existing allegedly LO, as well as NLO analyses of 
F2{x,Q'^) are incomplete. The source of this incompleteness of the conventional approach is 
traced back to the lack of clear identification of QCD effects and to the misinterpretation of 
the behaviour of q^{x,M) as a function of as(M). Complete LO and NLO QCD analyses of 
F2{x, Q^) are shown to differ substantially from the conventional ones. Whereas complete NLO 
analysis requires the knowledge of two so far uncalculated quantities, a complete LO one is 
currently possible, but compared to the conventional formulation requires the inclusion of four 
known, but in the existing LO analyses unused quantities. The arguments recently advanced 
in favour of the conventional approach are analyzed and shown to contain a serious flaw. If 
corrected, they actually lend support to my claim. 
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1. Introduction 



Observed from a large distance the photon behaves as a neutral structureless object governed 
by the laws of Quantum Electrodynamics. However, when probed at short distances it exhibits 
some properties characteristic of hadrons ^. This "photon structure" is quantified, similarly 
as in the case of hadrons, in terms of parton distribution functions (PDF), satisfying certain 
evolution equations. Because of the direct coupling of photons to quark-antiquark pairs these 
evolution equations are, contrary to the case of hadrons, inhomogeneous. This inhomogeneity 
has important implications for QCD analysis of -F^(a;, Q^) and other physical quantities involving 
photon in the initial state. 

In the previous paper [4] these implications led me to the conclusion the all existing NLO 
QCD analyses of F2{x,Q'^) are incomplete. With the exception of the authors of the FKP 
approach [5] my claim has been either ignored of rejected. The inertia of the "common wisdom" is 
enormous. I have therefore appreciated the recent attempt of A. Vogt [6] to present mathematical 
arguments in favour of the conventional approach. I analyze them in Section 4 and show that 
they contain a serious flaw. If corrected, Vogt's arguments actually lend support to my claim. 

In the course of discussions of this and related subjects concerning conventional QCD anal- 
ysis of photon structure, I have realized that the main source of the confusion surrounding the 
QCD analysis of -F2^(x, Q^) was the lack of a clear definition of what "LO" and "NLO" means in 
the context of photonic interactions. The point is that all existing conventional QCD analyses 
of F2 mix purely QED effects, which are usually quite dominant, with the genuine QCD ones, 
which represent mostly small corrections only. I will therefore start by presenting my definition 
of what "leading" and "next -to-leading" order means for parton distribution functions of the 
photon and for F2{x,Q'^). In the next step I will construct explicit solutions of the inhomoge- 
neous evolution equation including inhomogeneous as well as homogeneous splitting functions up 
to order a^. Straightforward analysis of these solutions shows that my claim in [4] was actually 
only partially correct: not only the existing supposedly NLO analyses are incomplete, but so are 
the LO ones! Contrary to the NLO analysis, which is currently impossible to perform because 
the necessary quantities have not yet been calculated, complete LO QCD analysis of F2{x,Q^) 
is feasible, but requires the inclusion of several additional 0{as) quantities. 

The paper is organized as follows. In the next Section basic facts and notation concerning 
PDF of the photon are recalled, followed in Section 3 by the discussion of the properties of 
the pointlike part of quark distribution function. Section 4 contains critical analysis of Vogt's 
arguments in [6]. In Section 5 the leading and next-to-leading order QCD analysis of F2 is 
formulated and the pointlike solutions of the inhomogeneous evolution equation for q^g are 
explicitly written down up to order q^. Phcnomcnological implications of the present analysis 
for F2 {x,Q'^) are outlined in Section 6, followed by the summary and conclusions in Section 7. 
Compared with [4] I have omitted the analysis of the structure of the virtual photon, some of 
which can be found in [7, 8]. 



^For recent theoretical and experimental reviews see [1] and [2, 3], respectively. 
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2. Notation and basic facts 



In QCD the coupling of quarks and gluons is characterized by the renormahzed colour cou- 
pling ( "couplant" for short) as , depending on the renormalization scale fi and satisfying the 
equation 

dW? = /3(a.(/^)) = -^aM - + • • • , (2-1) 

where, in QCD with Uf massless quark flavours, the first two coefficients, /?o = 11 — 2n//3 
and (3i = 102 — 38n//3, arc unique, while all the higher order ones are ambiguous. As we 
shall stay in this paper within the NLO QCD, only the first two, unique, terms in (2.1) will 
be taken into account in the following. However, even for a given r.h.s. of (2.1) its solution 
as{fi) is not a unique function of /U, because there is an infinite number of solutions of (2.1), 
differing by the initial condition. This so called renormalization scheme (RS) ambiguity ^ can 
be parameterized in a number of ways. One of them makes use of the fact that in the process 
of renormalization another dimensional parameter, denoted usually A, inevitably appears in the 
theory. This parameter depends on the RS and at the NLO even fully specifies it. For instance, 
as{n) in the familiar MS and MS RS are two solutions of the same equation (2.1), associated 
with different Ars ^- In this paper we shall work in the standard MS RS of the couplant. 

"Dressed" PDF ^ result from the resummation of multiple parton collinear emission off 
the corresponding "bare" parton distributions. As a result of this resummation PDF acquire 
dependence on the factorization scale M. This scale defines the upper limit on some measure t 
of the off-shellness of partons included in the definition of D(x, M) 

Diix,M)= / dtdi{x,t), i = q,q,G, (2.2) 

^min 

where the unintegrated PDF di{x,t) describe distribution functions of partons with the mo- 
mentum fraction x and fixed off-shellness t. Parton virtuality r =| — \ or transverse 
mass = + m?, are two standard choices of such a measure. Because at small t, 
di{x,t) = 0{l/t^),k = 1,2, the dominant part of the integral (2.2) comes from the region 
of small off-shellness t. Varying the upper bound in (2.2) has therefore only a small ef- 
fect on the integral (2.2), leading to weak (at most logarithmic) scaling violations. Factorization 
scale dependence of PDF of the photon ^ is determined by the system of coupled inhomogeneous 
evolution equations 



dS(x,M) 
dlnM2 
dG{x, M) 
dlnM2 
dgNs(x,M) 
dlnM2 



= 5^kq + Pqq®^ + PgG G, (2.3) 

= kG + PGq®^ + Pgg ® G, (2.4) 

= (^NS^ + ^*NS ® (2.5) 
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In higher orders this ambiguity includes also the arbitrariness of the coefficients l3i,i >2. 
^At the NLO the variation of both the renormalization scale /i and the renormalization scheme RS={Ars} is 
redundant. It suffices to fix one of them and vary the other, but I stick to the common practice of considering 
both of them as free parameters. 

*In the following the adjective "dressed" as well as the superscript "7" will be dropped. 
^If not stated otherwise all distribution functions in the following concern the photon. 
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where ^ 



1=1 i=l 
q^six,M) = ^ (e^ - (e^)) {qi{x, M) +qi{x, M)) , 



i=l 

To order a the sphtting functions Pij and fcj are given as power expansions in as{M): 

2 



kq{x,M) = 



a 
2^ 



a 



kG{x,M) = — 



Pij{x,M) 



27r 



X +••• 



27r 

„(0) 



27r 



(2.6) 
(2.7) 
(2.8) 

(2.9) 
(2.10) 
(2.11) 



where the leading order sphtting functions k^\x) = {x^ + (1 — xY) and P-j'{x) are unique, 



(0), 



^3 



while all higher order ones kq'\ k^Q, Pjfi\j ^ 1 depend on the choice of the factorization scheme 
(FS). The equations (2.3-2.5) can be recast into evolution equations for qi{x , M) ,q,i{x , M) and 
G{x,M) with inhomegenous splitting functions k^^ = 3ejkq^\ The photon structure function 
F2{x,Q^), measured in deep inelastic scattering of electrons on photons is given as 



a 



-F^{x, Q') = q^s{M) ® CgiQ/M) + — <5nsC^ + 

X ZTT 

a 



(2.12) 



e")S(M) ® Cq{Q/M) + ^{e')5^C^ + {e^)^G{M) ® Cg{Q/M) (2.13) 



27r 



27r 



of photonic PDF and coefficient functions Cq{x),CG{x),C-y{x) admitting perturbative expan- 
sions 



C,(x, Q/M) = 5(1 -x) + ^^CfHx, Q/M) + • • • , 



27r 

^C«(x,Q/M) + ..., 
C^{x, Q/M) = Cf (x, Q/M) + ^C(i)(^, Q/M) + ■■■, 



Cg{x,Q/M) = 



(2.14) 
(2.15) 
(2.16) 



where the standard formula for c!^^ reads ^ 



C^o)(x,Q/M) = (x2 + (l-x)2) 



+ 8x(l -x) - 1. 



(2.17) 



For nonsinglet quark distribution function qns another definition is used in the hterature [9]. The definition 

adopted here corresponds to that used in [10]. 

^Alternatives to this expression for = are discussed in Subsection 6.6. 
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The renormalization scale ^, used as argument of as(/u) in (2.14-2.16) is in principle indepen- 
dent of the factorization scale M. Note that despite the presence of /x as argument of as{n) 
in (2.14-2.16), the coefiicient functions Cq^Cc and C-y are actually independent of ^ because 
the //-dependence of the expansion parameter agifJ') is cancelled by explicit dependence of 

(i) (i) (i) 

Cq ,Cq ,C-y ,i > 2 on /X [11]. On the other hand, PDF and the coefficient functions Cq^Ca 
and do depend on both the factorization scale M and factorization scheme, but in such a 
correlated manner that physical quantities, like Fq^., are independent of both M and the FS, 
provided expansions (2.9-2.11) and (2.14-2.16) are taken to all orders in as{M) and q;s(//). In 
practical calculations based on truncated forms of (2.9-2.11) and (2.14-2.16) this invariance is, 
however, lost and the choice of both M and FS makes numerical difference even for physical 
quantities. The expressions for C'^q \ Cq^ given in [12] arc Tisually claimed to correspond to "MS 
factorization scheme". As argued in [13], this denomination is, however, incomplete. The adjec- 
tive "MS" concerns exclusively the choice of the RS of the couplant as and has nothing to do 
with the choice of the splitting functions P^p . The choices of the renormalization scheme of the 
couplant as and of the factorization scheme of PDF are two completely independent decisions, 
concerning two different and in general unrelated redefinition procedures. Both are necessary 
in order to specify uniquely the results of fixed order perturbativc calculations, but we may 
combine any choice of the RS of the couplant with any choice of the FS of PDF. The coefficient 
functions Cq.Cc.C-^, depend on both of them, whereas the splitting functions depend only on 
the latter. The results given in [12] correspond to MS RS of the couplant but to the ""minimal 
subtraction" FS of PDF ^. It is therefore more appropriate to call this full specification of the 
renormalization and factorization schemes as "MS -|- MS scheme" . Although the phenomeno- 
logical relevance of treating ^ and M as independent parameters has been demonstrated [15], I 
shall follow the usual practice and set fi = M. 

3. Pointlike solutions and their properties 

The general solution of the evolution equations (2.3-2.5) can be written as the sum of a partic- 
ular solution of the full inhomogeneous equation and the general solution of the corresponding 
homogeneous one, called hadronic ^ part. A subset of the solutions of full evolution equations 
resulting from the resummation of series of diagrams like those in Fig. (1), which start with the 
pointlike purely QED vertex 7 — ^ qq, are called pointlike (PL) solutions. In writing down the 
expression for the resummation of diagrams in Fig. 1 there is a freedom in specifying some sort 
of boundary condition. It is common to work within a subset of pointlike solutions specified by 
the value of the scale Mq at which they vanish. In general, we can thus write {D = q, q, G) 

D{x, M^) = D^^{x, M^) + D^^^{x, M^). (3.1) 

Due to the fact that there is an infinite number of pointlike solutions q^^{x, M^), the separation 
of quark and gluon distribution functions into their pointlike and hadronic parts is, however, 
ambiguous and therefore these concepts have separately no physical meaning. In [7] we discussed 
numerical aspects of this ambiguity for the Schuler-Sjostrand sets of parameterizations [16]. 

*Scc Section 2.6 of [14], in particular eq. (2.31), for discussion of this point. 

''Sometimes also called "VDM part" because it is usually modelled by PDF of vector mesons. 
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Figure 1: Diagrams defining the pointlike parts of NS quarlc distribution function in LL approximation. 
The resummation involves integration over parton virtualities Mg < r < M^.. 



To see the most important feature of the pointhke part of quark distribution functions that 
will be crucial for the following analysis, let us consider in detail the case of nonsinglet quark 
distribution function q^six,M), which is explicitly defined via the series 



27r Jm2 t Jx y \yj J A 

Jx y 



dri asin) a (o). . dr2 



p(0) ± 



w '''' \wy Jm^ 



dTiasin) r' dT2as{T2) a (o) 



Mi 



IM§ Ti 27r Jm§ T2 27r 27r 
where A;^g(x) = SNskg'\x). In terms of moments defined as 

/(n) = /' x^f{x)dx 
Jo 



■NS 



{w) / 



+ 



(3.2) 
(3.3) 



this series can be resummed in a closed form 

in 



q^^{n,Mo,M) 



as{M) 



where 



I 



a 



a, (Mo) 



l-2P,(^'(n)//3o 



aNs("'), 



(3.4) 



(3.5) 



27r/5oi_2P,(°)(n)//3o' 

It is straightforward to show that (3.2) or, equivalently, (3.4), satisfy the evolution equation 
(2.5) with the splitting functions kg and Pij including the first terms fcg^^ and Pqq^ only. 

Transforming (3.5) to the x-space by means of inverse Mellin transformation we get aNs(2;) 
shown in Fig. 2. The resummation softens the x— dependence of aNs(a;) with respect to the 
first term in (3.2), proportional to k-Ms{x), but does not change the logarithmic dependence of 
^NS on M. In the nonsinglet channel the effects of gluon radiation on g^g are significant for 
X > 0.6 but small elsewhere, whereas in the singlet channel such effects are marked also for 
X < 0.5, where they lead to a steep rise of xq^^ at very small x. As emphasized long time ago 
by authors of [5] the logarithmic dependence of q^g on InM has nothing to do with QCD and 
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Figure 2: Comparison, left in the nonsinglet and right in the singlet channels, of the functions xki{x) 
x5^kf\i = NS, S with the funct 
analogue in the singlet channel. 



xSik^\i = NS,S with the functions xa{x) corresponding to asymptotic pointlike solutions (3.6) and its 



results exclusively from integration over the transverse momenta (virtualities) of quarks coming 
from the basic QED "f ^ qq splitting. For M/Mq oo the second term in brackets of (3.4) 
vanishes and therefore all pointlike solutions share the same large M behaviour 

q^^ix, Mo, M) ^ -^a^sioc) = qMi^, M) oc In (3.6) 

defining the so called asymptotic pointlike solution q^l{x,M) [17, 18]. Note that (3.6) can be 
interpreted as a special case of the general pointlikc solution (3.4) in which the lower integration 
limit Mq in (3.2) is set equal to A, i.e. Mq = A! The fact that for the asymptotic pointlike 
solution (3.6) as{M) appears in the denominator of (3.6) has been the source of misleading 
claims (see, for instance, [1]) that q(x,M) = 0{l/as). This claim is wrong for a number of 
reasons. First, it is manifestly invalid for the widely used Schuler-Sjostrand SaSl and SaS2 sets, 
which take Mq = 0.6 GeV and Mq = 2 GeV. It is obvious [19] that provided Mq is kept fixed 
when as ^ the sum (3.2) approaches its first term, i.e. 

q^^ix,M,Mo) ^ ^4°)(x)ln^, (3.7) 

corresponding to purely QED splitting 7 — ^ qq. However, the claim that g oc l/a^ is misleading 
even for the asymptotic pointlikc solution (3.6). The fact that it diverges when A is a direct 

consequence of the fact that for this (and only this) pointlike solution the decrease of the coupling 
as{M/A) as A ^ is compensated by the simultaneous extension of the integration region as 
Mq = a — > 0! If, however, QCD is switched off by sending A — > without simultaneously 
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extending the integration region, i.e. for fixed AIq, there is no trace of QCD left and we get 
back the simple QED formula (3.7). This suggests to identify QCD contribution to NS quark 
distribution function with the difference 



ggg^^(x, M) ^ q^^{x, Mo,M) - ^k'Slix) In ^. (3.8) 



^A;(°VxUn — 



Expanding (3.4) in powers of as{M), keeping Mq and M fixed we find 

as expected from the explicit expression for the second term on the r.h.s. of (3.2). Note that the 
leading term in (3.9) is proportional to the second power of Lm = ln(M^/MQ). In subsection 
5.2.2 I will argue that in a complete LO QCD analysis qns contains another term that behaves 
as agLM- 

In summary, q^^{x,M) can be written a sum of two terms: purely QED contribution, 
proportional to a, and the term describing genuine QCD effects, proportional to aag and thus 
vanishing when QCD is switched off. One can use the usual claim that q oc a/ag merely as a 
shorthand for the specification of large M behaviour of (3.2), as expressed in (3.6). In fact, this 
is in fact what one finds in the original papers [17, 18]^° which do not contain any explicit claim 
that q{x,M) = 0(1/ ag). 

4. What is proving Vogt's proof? 

Before presenting my reformulation of the LO and NLO QCD analysis of F2 , let me go through 
Vogt's arguments that purport to prove that q{M) oc a/ag. For the purpose of the discussion 
in this Section I will adopt Vogt's simplified notation, in which Ug = ag/Air, F'^ = F^-^^/a, (p = 
Ins dF^/dlnQ'^ = , all obvious indices as well as overall charge factors are suppressed 
and products of x-dcpcndcnt quantities are understood as convolutions in x space or simple 
products in momentum space. We have 

F^ = Cgq^ + C^ (4.1) 
f = k + Pq^ (4.2) 
FT = Cqq^ + Cgf + (4.3) 

The r.h.s. of (4.3) can be rewritten as [6] 

F^ = (^-CgCg ^C^ + Cgk - PC^ + C^) + [CgC-^ + p) F\ (4.4) 

where expressions in the brackets are factorization scheme invariants. Inserting into (4.4) per- 
turbation expansions of Cq,C~^, their derivatives and the splitting functions k, P we get 



i°See, for instance, eq. (3.20) of [18]. 

^^Note that due to slightly different notation my Cif^ plays the role of Cj^'^ in [6]. 



^^Note that eq. (14) of [6] contains two misprints in the bracket standing by a^: the last term comes with the 

)ng sign and the 
in Vogt's notation). 



wrong sign and the preceding one should correctly read —fioC^^ (my C^^^ and C^^ correspond to C^^^ and 
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+{asf [fc(^) + C(^)A;W - pWcf ] (4.5) 
+FTKP(0) + •••]. 

The fact that C^^ appears in the expression by Og in (4.5) together with in Vogt's 

words that "C^'^^ .... enters F'^ on the same level as the hadronic NLO quantity Cg^^" leads him 
to the conclusion that "C^°^ is to be considered as a NLO contribution", and, consequently, " 
for the purpose of power counting in ag the quark densities and the structure function F2 
have to counted as l/a^." The flaw of this argument is obvious if one applies it to C^^^ (^7^'' 
in Vogt's notation) appearing in the expression standing in (4.5) by Og accompanied by both 
and C''q \ Repeating Vogt's argument leads us to contradictory conclusions that C^"^ is 
simultaneously of the same order as C^^^ and ! The resolution of this contradiction is obvious: 
we have to take into account the fact that C^°^ and C-y^^ do not enter the coefficients in (4.5) 
alone, but in products with other quantities, corresponding to different orders of perturbative 
QCD. For instance, since k^^^ stands in (2.9) by {asf* whereas k^^^ by ctg, the products C^'^k^^'^ 
and Cq arc of the same order despite the fact that c'^'^ and c'^q^ stand in (2.14) by 

as and respectively. Similarly, C^^ is, as expected, of one order of higher than C'^\ 
The failure of Vogt to take this fact into account led him to wrong conclusion 
Taking into account that k^^^ stands by {asT whereas P^^^ by ag implies that Cq^^ is of one 
order of ag higher than C-y^^ and thus C^^ not of the NLO! 

If I am right, why have most theorists so tenaciously clung to the claim that in some 
sense q cc l/ugl In part because it provides a simple way of expressing large scale behaviour 
of q{x,M), but the main reason is related to another tenet of the conventional approach to 
F2{x,Q'^), namely the assumption that at the leading order of QCD the observable F2{x,Q'^) 
is related to q{x, M) in the same way as for hadrons, i.e. 

Fl^0^x,Q^) = qi^o{x,M^). (4.6) 

This - incorrect as I am just going to argue - relation is in fact the principal source of all 
misleading and factually wrong statements concerning the QCD analysis of F2. 

To evaluate the r.h.s. of (4.6) the factorization scale M has to be chosen. The requirement of 
factorization scale invariance of physical quantities implies that at any finite order of perturbative 
QCD the difference PQ^"'\Mi) — PQ^^\M2) of perturbative predictions for a physical quantity 
PQ evaluated up to order a" at two scales Mi and M2 behaves as 0(q;"^^). For hadrons (4.6) 
is consistent with this fundamental requirement due to the fact that for them q{x, M) satisfies 
homogeneous evolution equations. Consequently, the difference 

A(Mi,M2) = q{Mi)-q{M2) oc a,(Mi)ln (^^J g(Mi) (4.7) 



Which is one order of Os higher than Cq^^kq°K 
^''With exception of the authors of [5], who have advocated ideas closely related to those advanced in this paper, 
for a long time. 
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is, indeed, of one order of a<j higher that q itself and therefore (4.6) to the order considered 
independent of the choice of M\ 

For the photon the presence of the inhomogeneous term (a/27r)feq°^ in the evolution equations 
for quark distribution functions implies, however, A(Mi, M2) cx a In (Mf/Mf) and, consequently, 
F^lq(Mi) — F^-^q{M-\) (X a. To retain the factorization scale invariance of (4.6) for -F^lq 

property q oc l/cts seems therefore indispensable! The use of the term "NLO" for C^°^ is also 
sometimes justified by the fact that for fixed ratio d = M1/M2, the dependence of difference 
A{Mi, M2) oc alnd is "subleading" with respect to gLO oc alnM^. This is true but it must be 
kept in mind that the dominant In behaviour of the standard definition of ^lo is basically 
a consequence of QED, not QCD dynamics! To claim with Vogt that the purely QED term 
Cj'^ (Q^ = M^) is of the "NLO" merely because it is M^-independent constant and thus small 
for large with respect to (/qed oc a In is a misuse of the terms "LO" and "NLO" , which 
are meant to describe different orders in a^, not different large behaviours. 

In the conventional approach to F2 wrong medicine is thus used to salvage the factorization 
scale invariance of (4.6). There is, however, a different and consistent way of guaranteeing this 
invariance at the LO that does not rely on the untenable assumption g oc l/a^: to modify the 
relation (4.6) itself! The rest of this paper is devoted to the elaboration of this claim. I will 
show that the modified relation between I^i^q and gLO is consistent with the obvious fact that 
g oc a and demonstrate how the additional terms in this relation conspire to quarantee the 
factorization scale invariance of F^^q. 

5. QCD analysis of F^(x, Q^) 

In this Section QCD analysis of F2{x^Q'^) that separates genuine QCD effects from those of 
purely QED origin will be presented. Throughout this and following sections I will restrict my 
attention to the nonsiglet part (2.12) of F^ 

Q^) = qm{M) ® Cg{Q/M) + ^S^sC^ (5.1) 

and nonsinglet quark distribution function as defined in (2.7). I will consider the contribu- 
tion of light quarks (i.e. u, d, s) only and, moreover, disregard the difference between their 
distributions functions after division by e^. Under these simplifying assumptions we have 

-F^^^si^, Q") = (5ns \q{M) Cg{Q/M) + ^C^{Q/M)\ (5.2) 

= 5ns [qiM) + ^g(M) ^ C«(g/M) + ^cf (Q/M) + ^^c\^)(Q/M) ■ ■ ■ 

where q = u{x)/3e'^ = d{x)/3e'^ = s{x)/3e^. For simplicity I will in the following drop the 

overall charge factor (5ns SiS well as the subscript "NS" and continue to use the dot to denote 

the derivatives with respect to InM^. In (5.1-5.2) I have written out explicitly the symbol 

to denote convolutions in x-space, but henceforth I will work mostly in the momentum space 

and thus all products of quark distribution and coefHcient functions are understood as simple 

multiplications in momentum space 

^^This difference is entirely negligible above roughly x = 0.25 but becomes sizable below this value. 
^®To save space, their dependence on the momentum variable n will not be written out explicitly. 
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5.1 Defining leading and next— to— leading orders for F2(x, Q^) 

Let us first clearly define what is meant under the terms "leading" and "next-to-leading" orders 
in the case of QCD analysis of As for hadronic parts of photonic PDF the definition of these 
terms is the same as for hadrons, I will concentrate in this Section on the properties of the 
pointlike quark distribution function and its contribution to 

It is useful to recall the meaning of the terms "leading" and "next-to-leading" order of 
perturbative QCD for the case of the familiar ratio 

measured in e'''e~ annihilations. The prefactor 

i?QED^ (^S^^e?^ (5.4) 

comes from pure QED, whereas genuine QCD effects give r{Q) as perturbation expansion in 
powers of ag 

as{M) 

TT 

For the quantity (5.3) it is a generally accepted procedure to divide out the QED contribution 
-Rqed and apply the terms "leading" and "next-to-leading" only to QCD analysis of r{Q), which 
starts as ag/n. Nobody would suggest including -Rqed in the definition of the term "leading" 
order in QCD analysis of (5.3). 

Unfortunately, precisely this is conventionally done in the case of F^. I will now present the 
organization of QCD expression for F2 that follows as closely as possible the convention used 
in QCD analysis of (5.3). For this purpose let us write, following the discussion at the end of 
Section 3, the pointlike quark distribution function q{M), satisfying the condition q{Mo) = 0, 
as a sum of two terms 

q{M) = (Zqed(M) + gQCD(M) (5.6) 
where the purely QED contribution is given as 



riQ) 



l + ^n(Q/M) + 



(5.5) 



9QED(M)^^fef In^. (5.7) 

Mq, a free parameter specifying the pointlike solution, can be interpreted as the lower limit on 
the integral over quark virtuality included in the definition of qQEDiM). Defined in this way 
9QCD is due entirely to QCD effects, i.e. vanishes when we switch off QCD. The expression 
(5.7) is a close analogue of the QED contribution .Rqed in (5.4). Note that qqcd satisfies the 
inhomogeneous evolution equation 

£^S^^QCD + (^)Vi^QCD + ---, (5.8) 
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which differs from that satisfied by the full quark distribution function not only by the absence 
of the first term {a/2Tr)kq^^ but also by shifted appearance of higher order coefficients kq^; i > 1. 
For instance, the inhomogeneous splitting function A;^^^ enters the evolution equation for qqcd 
at the same order as homogeneous splitting function Pgq^ and thus these splitting functions 
will appear at the same order also in its solutions. The fact that kq^^ is a function of x (or 
n) only, whereas qqeb is in addition also a function of M, influences relative importance of 
the two terms making up the coefficient by in (5.8) as M varies, but docs not change the 
basic observation, namely that both of them contribute at the same order. Similar statement 
holds for all pairs kq~^^\Pgg;i > 0. The simultaneous presence of k^^ and Pql^ in the O(a^) 
term of the inhomogeneous part of (5.8) is yet another expression of my claim in [4] that NLO 
QCD analysis of F2 requires the knowledge of k^\ I have, however, failed to realize that 
the same argument implies that in a complete LO QCD analysis the evolution equation for 
quark distribution function must include the 0{as) inhomogeneous splitting function kq^^ as 
well. Explicit expressions of the resulting solutions are presented and their properties discussed 
below. 



5.1.1 QED part 

The purely QED result for F2 is given by the sum 

^^2^QED(a^> Q^) = «Qed(x, M) + ^C(o)(x, Q/M), (5.9) 

which is a close analogue of QED prediction (5.4) for the ratio (5.3). Note that although both 
terms on the r.h.s. of (5.9) depend on factorization scale M, their sum is independent of it and 
is a function of and Mq only! 

Let me emphasize that already at this point I depart from the conventional approach which 
identifies df^ as part of the "NLO" corrections. In fact df^ has nothing to do with QCD 
at all and is entirely due to QED effects! It must therefore be always present in any QCD 
analysis of data. It is thus also not true, as claimed in [20], that the SaSlM and SaS2M 
parameterizations arc "theoretically inconsistent" because they combine in LO expression 
for F2 the "NLO" quantity c!f^ with the LO quark distribution functions. In fact, just the 
opposite is true! Prom the point of view of including dy*^ the analysis of F2 in terms of SaSlM 
and SaS2M parameterizations are closer to complete LO QCD analysis than analogous analysis 
using SaSlD or SaS2D sets. 

Let me reiterate that the principal feature of F2{x,Q'^), namely its logarithmic rise with 
, is basically a QED effect and consequently its observation in experiments by itself no check 
of QCD. One has to go to subtler features '^^ to identify genuine QCD effects! 

'^'^In which the integration over virtual quark virtualitios is restricted to r > Mq . 
^*The dependence on Mo must not be confused with factorization scale dependence. 

^®In fact the SaSlM and SaS2M parameterizations include in only part of the expression (2.17). This and 

related facts are discussed in detail in Subsection 6.6. 

^°For instance, the deviation of the x-dependence of dF^ {x,Q^)/d\nQ^ from the QED result {a/2-K)kq'^\x). 
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5.1.2 QCD part 

In conventional analyses of F2 the first term on the r.h.s. of (5.9) is included in the "leading" 
order of QCD thereby obscuring meaning of the term "leading" . There is no obstacle to following 
the procedure adopted in QCD analyses of (5.3) and rewrite also F^{x,Q'^) as the sum of its 
QED and QCD parts 

-i^2'' = ?QED + ^Cf+ (5.10) 
pure QED 

.QOO + gc«,«.o + ^g<7«+ (5.11) 

V ' 

=Ai, starting as 0{aas) 

£c...,,co.^(^)^Cf.(£)^C.^.«.o ,5.12) 



=j42, starting as 0{aaj) 

In (5.10-5.12) I have grouped into quantities Aq, A\ and A2 the contributions that start at zero, 
first and second order of cig, assuming gqcD ~ Q!s- In the next subsection I will present explicit 
solutions demonstrating this behaviour. 

The first line in (5.10) contains the purely QED contribution (5.10) to F^ - Following the 
analogy with (5.3) the LO QCD approximation to F^ is identified with A\, 

-/l^oM-) = ,QCD + ^C^'W^ + (5.13) 

the NLO one with the sum Ai + A2, etc.. 

The expression (5.13) represents an analogue of the LO contribution to (5.5), equal to 
''LO = cts/n. Beside its more complicated structure eq. (5.13) differs from tlo also by the fact 
that whereas the latter is unique, (5.13) depends on Mq. This means that any QCD analysis of 
F2 must start with fixing the value of this parameter. It is worth emphasizing that for heavy 
quarks Mq can be related to uiq and for virtual photons to their virtuality P^. 

As my criticism in Section 4 of the conventional formulation of QCD analysis of F^ was 
based on the lack of the factorization scale invariance of the latter, let us now check whether 
this invariance holds for the expression (5.13) where ^qcd satisfies (5.8) up to order a^. Taking 
the derivative of (5.13) with respect to we find, taking into account that from similar 
considerations in hadronic collisions we know that Cg^^ = —P^ , 



Tpi _ a as 
^LO-27r27r 



+ 0{al). (5.14) 



Contrary to the case of hadronic structure function Fj', where analogous derivative is propor- 
tional to asF2 and thus manifestly of NLO, for (5.13) the condition that (5.14) is of order 
implies the following nontrivial relation 



4'' + CJ'^ + = (5.15) 
'j!y^\ Because 

actually implies two nontrivial relations. Their validity can be verfied using the expressions for 



between the quantities kq^\ kq^^ and C^y^\ Because c!^^ depends on InM^ as ln^(QV^^) (5-15) 



12 



NLO coefficients Cg^^ and Cq \ calculated in [21, 22]. According to [23] and taking into account 
slightly different normalization convention C^^^ can be obtained from Cq^ in [22] by 

• replacing TfUf by unity, 

• dropping terms including Pqq and Po, which are absent for Cj^\ and 

• replacing P^^ with kq^^ . 

Substituting c!^^ obtained in this way into (5.15) one finds that the sum on its l.h.s. indeed 
vanishes. Let me emphasize that for this cancellation the presence of kq^^ in (5.15) is vital. 

For the purpose of comparing LO QCD expression for F2 , defined in (5.13), with the 
conventional LO formula as well as with the data, we add to it the QED contribution (5.9): 



X V 



^QED + ^2%) = 5QED + ^QCD + ^Cf^ + gc(l),QED + ^^C^'^ ■ (5-16) 



Recalling the conventional formulation of the LO QCD approximation of F2 [6] 

iF2%(x,Q2) = g(x,Q), (5.17) 

where q(x, M) satisfies the evolution equation that includes kq'^ and splitting functions 
only, we see that it differs from (5.16) 

• by the absence of the contributions of photonic coefficient functions C^^ and C^^^ , 

• by the absence of the convolution of quark coefficient function C^^ with qqed, and 

• by the fact kq'^ is not included in the evolution equation for q{M). 

All these differences are substantial, but as all necessary quantities are available, there is no 
obstacle to performing complete LO QCD analysis using formula (5.16) with gqcD given by the 
formula (5.32) of the next subsection. 

On the other hand, as neither kq ' nor are known, a complete NLO QCD analysis of 
P^ is currently impossible to perform. Within the conventional approach the terms in (5.13) 
proportional to and are part of the NLO expression 

^^2Vo = ^+i^i^^^ + £c(°), (5.18) 

but the photonic coefficient function though known, is not used even at the NLO. As we 
shall see in Section 6, its contribution is numerically comparable to that of C^^\ 



^^My are by a factor of 4 smaller than those used in [21, 22], and Pij are expanded in powers of a/27r in 
this paper and in powers of Qs/47r in [21, 22, 23]. 
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5.2 Explicit form of the pointlike solution qNs(n, M, Mq) 

The inhomogeneous evolution equation for nonsinglet quark distribution function is technically 
particularly simple to solve in the case (3i = 0,i > 1. As none of the conclusions of this paper 
depends in any essential way on nonzero value of Pi and higher order coefficients > 2 
can be set to zero by the choice of an appropriate renormalization scheme, 1 will set (3i = 
0;i > 1 thorought the rest of this paper. Under this simplifying assumption dg(M)/dlnM^ = 
(— /3o/47r)dg/das and the inhomegeneous evolution equation can be rewritten as 



I'iy) ^ 



dq{y) _ 47r 1 dq{M) _ 
dy ~ ~A^afdlnM2 " 

1 aki'^ 1 
+ —7^ h 



2ak, 



(3o y2 7r/3o y 27r2/?o 
-Qiy) - Piy)qiy) 



+ 



9p(0) 1 p 



(1) 



-- + 



+ 



Pr 



(2) 



Po y 7r/3o 27r2/3o 



-y + 



where 



9p(0) 1 p{l) 



+ 



P 



(2) 



f3o V vr/3o 27r2/3o 



-y 



i=l 



Qiy) 



2aA:(0) 1 ak^^^^ 1 aA:(2) 

+ — + TT^TTT + 



Po y^ nf3o y 27r2/?o 
This standard inhomogeneous differential equation of the type 

q'{y) + P{y)Q{y) + Q{y) = o 

has a general solution in the form 



i=l 



Qiy 



i-3 




(5.19) 

(5.20) 
(5.21) 

(5.22) 




Figure 3: The functions V{z,p) (solid) and W{z,p) (dotted), where p = —2Pqq\n)/po, evaluated for 
n = 1, 3, 6 and conipared to f{z) = z (dashed lines). 



. X - P(u)du \ r r^f \ f P(u)du 

q{y,yo) = e ■'yo ^ ' go- Q{v)&yo ^' 

Jyo 



(5.23) 
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where qq = q{yo) = qiMo) specifics the initial condition imposed on the solution of (5.23). Our 
pointlike solutions correspond to go = and we thus have 



iiy,yo) 



- r Piu)du 



r Q{v)e^yo = - r Q{v)e^: ^W^", (5.24) 



'yo •Jyo 
For further discussion it is useful to introduce the following functions 



W{z,p) = 



P 
-2 



p(l +p)z 



[1 - (1 - z^+P - (1 + p)i 



(5.25) 
(5.26) 



where 



z = 



y yo 



=yA = as{M)^LM = 



In 



M2 



In 



M2 



Lm = In 



(5.27) 



which both behave for z ^ and all p like f(z) = z. Their shapes are shown, for three values 
of p corresponding to n = 1,3,6, in Fig. 3. For fixed M^,Mq the region z ^ corresponds to 
the limit ag = y ^ 0. Inserting into (5.24) the expansions (5.20) and (5.21) and retaining in 
both cases first three terms we get 

q{y,yo) = e-(P'y+-^P'y") 

q,Viz,l-pi)- + q,p^V{z,-pi) + SM±pAv{z, -1 - pi)y 

y ^ 



+ q2V{z,-pi) + q2P2V{z,-l - pi)y + 



q2{pl +P3) 



y(z,-2-pi)2/2 



+ q3Viz,-l-pi)y + q3P2Viz,-2-pi)y^ + ^^^^^^Viz,-3-pi)y^ 



(5.28) 



Expanding further the exponential prefactor in (5.28), adding and subtracting the purely QED 
contribution (5.7) and grouping together term that stand by same power of y = as{M) we get, 
keeping terms up to y^ = a^, 



+ 



+ 



ly 



V{z, 1 - pi) - — In ^ + q^Viz, -pi] 



Q1P2 



V{z, -p^)-V{z,l-pi) 



y 



+ qsViz, -1 -pi) - q2P2V{z, -pi) 



(5.29) 
(5.30) 
(5.31) 



The first line in the above expression defines the purely QED part, Q'qeD) the second line the 
LO approximation, ^lO) a-^d the sum of (5.30) and (5.31) ^nlo- Note that the difference 
{V{—pi) — V{1 —pi))/y behaves at small y as \A^y = 0{y), i.e. in the same way as V{z,p) « 
z = yA = 0{y) and therefore contributes to ^nlo- 
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5.2.1 QED part 

The logarithmic rise of F2{x,Q'^) with results from integration of transverse momenta (vir- 
tualities) of quarks/antiquarks produced in the primary QED vertex ^ ^ qq and has therefore 
nothing to do with QCD. The scale Mq, which from the point of view of mathematics defines 
the scale at which the initial conditions on (5.22) are imposed, has a clear physical meaning: 
it separates the region of transverse momenta (virtualities) where nonperturbative effects are 
dominant [5]. Note that although ^qed and separately depend on the factorization scale 
M, their sum -P2''qed '^oes not! The purely QED part (5.9) of F2 plays similar role as the purely 
QED prefactor 3^e? in (5.3). We can either subtract it from data or include it in theoretical 
analyses. The second option, i.e. considering the full sum q{x, M) in (5.6) as the basic quantity, 
is preferable because the QCD part of q{M) alone, i.e. qqcD, is negative In the following we 
shall analyze the basic features of -P^^qcd and gqcD but keep in mind that for proper physical 
interpretation the purely QED parts F^q^^ and qqsD have to be included as well. 

5.2.2 QCD part: leading order 

In terms of functions V{z,p) and W{z,p) the LO expression for gqcD reads 



9lo(m2 



a k, 



(0) 



27r Po 



p(0) jyj2 



a 2k^ 



(1) 



-Viz, -pi). 



(5.32) 



The first term in (5.32) coincides with the LO expression (3.4) of the conventional approach 
from which the purely QED part (5.7) was subtracted. The importance of the second term 
in (5.32), which is absent in the conventional approach, will be discussed in the next Section. 
The order of (5.32) is determined by its behaviour in the limit 0. In investigating this 

behaviour all external kinematical variables, as well as factorization scales must be kept fixed 
and the approach — realized by sending Aqcd — > 0. In our case this means keeping Q^, 
and Mq fixed and sending z ^ 0. In this limit both terms in (5.32) behave in the same way and 
we get 



QLo(M) 



qiPiA"^ 



+ q2A 



as{M) 



a as{M) 
2^ 27r 



^(0) p{0) 



11 



Lli + k^}^L 



M 



(5.33) 



Eq. (5.33) is closely analogous to standard LO expression for r{Q) in (5.5). It is easy to check 
that (5.33) satisfies (5.8) where only the first term on its r.h.s. is taken into account. 

The behaviour of qi^Q[x,M) at large scales M is determined by (5.32) in the limit z — >■ 1. 
Taking into account that V{l,p) = 1/p and W{l,p) = 2/(1 + p) we find 



qi.o{x,M) 



,2R 



(0) 



11 



1 



27r ^ /3o l-2PgVA 



a 



k 



(1) 







27r_p(0)' 



(5.34) 



The first term in (5.34) is dominant as M — 00, but as we shall see in the next Section the 
second one, absent in conventional LO analysis, is actually numerically more important in most 
of the currently accessible kinematical region. 

■^■^Thc fact that QCD correction to q is negative is nothing extraordinary. The same holds, for instance, for 
QCD corrections to Gross-Llewellyn-Smith sum rule. 
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5.2.3 QCD part: next to leading order 

The expression for qNLoiM"^), given by the sum of (5.30) and (5.31), behaves at small y as 



gNLo(M2) 



qiPi 



^2 



+ q2A 



as{M) + 



+ (93 - 92^2)^ 



(5.35) 



as{M) 
27r 



1.(0) p(0) 



+ 



«s(M) 
27r 



91.(1) p(l)' 



L 



M 



and involves in addition to quantities entering the LO formula (5.32), i.e. klj^\ k^^ and Pqq' also 

the 0{af) inhomogeneous splitting function and 0(a'^) homogeneous splitting function Pql^ ■ 

Whereas the latter is known, the former is not, thereby preventing the evaluation of (J'nlo- In 

addition to this obstacle, a complete NLO analysis of FJ requires also the knowledge of O(aa^) 

(2) 

photonic coefficient function C-y , which is also so far unknown. Consequently, a complete NLO 
QCD analysis of F2 is currently impossible to perform. 



,(0) 



6. Phenomenological implications 

In this Section I will demonstrate numerical importance of the contributions to F^^^q resulting 
from the inclusion of the quantities that are omitted in standard LO analysis: the inhomogeneous 
splitting function kq^^ , photonic coefficient functions Cy^^ and Cy'^ and quark coefficient function 
Cq^\ All quantitative considerations concern the LO QCD analysis of F2 in the nonsinglet 
channel only, assuming, for the sake of technical simplicity, /3i = 0. Any realistic analysis of 
experimental data will require inclusion of effects of nonzero value of /3i, as well as extension 
of the formalism to the singlet channel and incorporation of contributions of hadronic parts of 
PDF. The work on these problems is in progress. Nevertheless, the quantitative impact of some 
of the effects discussed in this Section on the LO QCD analysis of F2 is so large that we may 
expect them to affect significantly the complete LO analysis of F2 as well. 

6.1 The effects of 

In all consideration of this Section I will take for the standard expression (2.17), evaluated 
for = Q^. The presence of the lower bound Mg on quark virtuality complicates, however, 
the situation and I will therefore discuss this point in more detail in subsection 6.6. 

6.2 The effects of C^^^ 

For illustration of the numerical importance of the 0{as) photonic coefficient function Cy^\ its 
contribution to LO analysis of F2 in MS+MS scheme is compared in Fig. 4 to those of C^^ {x) , 
the QED formula (5.7), as well as other QCD contributions discussed in this Section. The 
comparison is performed for = 10,25,100 and 1000 GeV^. Contrary to the contribution 
of c!f\ C-^^^ enters the expression for F^ multiplied by as{Q), and therefore decreases with 
increasing until it becomes negligible with respect to the former as Q ^ 00. However, one 
would have to go to extremely large Q^, inaccessible in current experiments, for this dominance of 
Cf to b e a good approximation. Fig. 4 shows that in the whole currently accessible kincmatical 
region of and for x up to about 0.85, the term proportional to C^^^ provides numerically the 
most important 0{as) correction to QED formula (5.9). 
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0.004 I 

- (a/27T)C'°' 

- '(a/2n)On(QVM^) 
0.003 - (ay27T)C%,p(Q) 

" ' first term in (5.32) 

" second term in (5.32) 

0.002 - 

(a/2TT)(a3/27T)C';' 

0.001 - 



-0.001 



0.005 



0.004 - 



0.003 - 



0.002 - 



0.001 - 



-0.001 - 




0,2 0.4 0.6 0.8 



0.2 0.4 0.6 



Figure 4: Comparison of individual contributions to F2 as described in the text. 



6.3 The effects of C 



(1) 



Performing tfie convolution C^^"* ® Q'qeDj implied by (5.16), where 



/ 9 + 5x 1 + 

l-x 



Inx — 



3 1 + x^ 

4 [1 - x\. 



+ (l + a;2) 



ln(l — x) 
l-x 



"3 



(6.1) 



is straightforward and leads to the thin solid curves in Fig. 4. This contribution has qualitatively 
the same shape as that of k^q \ discussed below, and is vital in the region close to a; = 1. 
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6.4 The effects of k^^ 




1 10 1 10 1 10 

[GeV'l [GeV'l [GeV'l 



Figure 5: The LO expression (5.32), multiplied by 2/3o7r/Q:, as a function of for three moments 
n = 1,3,6, Mo = 1 GeV and A = 0.2 GeV. The solid line shows the full result, the dotted and dashed 
ones correspond to the first and second term in (5.32) plotted separately. 

-2 

X 10 




Figure 6: Illustration of the approach to the asymptotic pointlike expression (3.6) (a), and the compar- 
ison of three expressions for c!y^\x) discussed in the text (b). 

At asymptotically large values of the first term on the r.h.s. of (5.32), which is the only one 
included in standard LO analysis of F2 , dominates due to the presence of "large log" ln(M^/MQ ) 
but again one would have to go to very large M^, inaccessible in current experiments to 
see this dominance numerically. This is illustrated in Fig. 5, where qi,o{M), as defined in 
(5.32), is plotted as a function of for three moments n = 1, 3, 6 and together with separate 
contributions of the first and second terms in (5.32). We see that for < 100 GeV^ the effect 

^^For Mo = 1 GeV and currently accessible < 200 GeV^ the "large log" ln(M^/M^) < 5, hardly a large 
number. 
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of including the latter almost caneels the negative contribution of the term appearing in the 
conventional analysis. Moreover, converting (5.32) into the x-space by means of inverse Mellin 
transformation we find (see Fig. 4) that the negative contribution of the first term in (5.32) to 
moments of qLo{n, M) comes mostly from the region close to x = 1. Below x = 0.75 the genuine 
QCD effects described by this part of (5.32) are tiny and smaller than those described by the 
second term, proportional to the inhomogeneous splitting function A;^^^! The latter gives small 
negative contribution to qijo{x,M) up to a; = 0.7 and large positive one above that value. Figs. 
4 and 5 show convincingly that in the whole region of and in most of the range of x accessible 
experimentally, the second term in (5.32) is numerically more important than the one included 
in standard LO analyses. Only for x close to x = 1 are both terms in (5.32) comparable and 
more or less cancelling each other. The tiny effect of QCD corrections described by the first 
part of (5.32) reflects the fact that in the region x < 0.65 its contribution to scaling violations 
of the sum qqed + <?QCD , plotted in Fig. 6a, are very small. 



6.5 Summary of QCD contributions at the LO 

In Fig. 7 the sum of contributions to LO analysis of F2 coming from inclusion of C^"^\ C^^"* 
and k^"^ is compared to that given by the flrst term in (5.32), the purely QED contribution 
(5.9), the complete LO QCD contribution (5.13), as well as full QED+QCD expressions in my 
(5.16) and standard (5.17) approaches, in the latter including the contribution of C^^ as well. 
In all calculations = Q^. The difference between the solid and dash-dotted curves, which 
represent my and standard expressions for i^^LO' large and persists for all experimentally 
accessible values of Q^. This is caused by the fact that in the whole kinematical region of 
and for most of the region of x accessible in current experiments, the sum of the three 
additional contributions included in the LO analysis proposed in this paper, is much larger than 
the contribution of the term generated by ■ But even for x close to x = 1 these additional 
terms are important as they help cancel in part the negative contribution of C^^ . 



6.6 Dispensing with the DIS-y factorization scheme 

The quantitative comparison in Fig. 7 reveals another interesting difference: the term ln(l — x) 
in (2.17) is much less troubling within my formulation than within the standard one. The 
reason is simple: the contributions of additional terms compensate large part of the negative 
contribution of C^^ at large x and shifts the region where the complete LO expression turns 
negative much closer to x = 1. 

A simple explanation of the origins of this term is provided within the parton model deriva- 
tion [4] of the expression (2.17) in which the parallel singularity associated with the splitting 
7 — > is regularized by means of quark masses and the troubling term ln(l — x) comes from 
the lower bound on the quark virtuality r (see Fig. 1) in coUinear kinematics 

^coU = '^l <T < -r^°ii - 

min — -j^ X 'max 

Performing the integration over r in the limits r^JJ < t < M'^ , as prescribed by the definition 



rr = 3^ < r < = ^- (6-2) 
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0,005 r 

; additional 0(aj terms 

conventional 0(a„) term 

0.004 - ^ 

: F,,,c.. eq. (5.13- 

- ' ^z.QED. eq. (5.9) 
0.002 ^ 

0,001 - - --**"'* 



: Q'=10GeV' 

-0.001 
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Figure 7: The sum of contributions to LO QCD expression for F2 included in this analysis, but omitted 
in the conventional ones, compared to the the effects of Pqq^ as well as to the QED formula (5.9), complete 
LO expression (5.16), and standard LO formula (5.17). 

(2.2) of PDF of the real photon, the primary QED splitting j ^ qq leads to [4] 

/(^) 1" {^) = /(^) 1- = /(^) 1- (^) + /(^) 1^(1 - (6-3) 

where f{x) = + (1 — x)^. Normally, the term proportional to ln(M^/mq) is included in the 
quark distribution function q{x, M), whereas the "constant" part of (6.3), i.e. f{x) ln(l — x) goes 
into the coefficient function c!f\ Note that we cannot set niq = in the logarithmic term. The 
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complete expression for nonlogarithmic term of (2.17) includes, beside the term /(x)ln(l — x), 
also terms coming from the upper bound in (6.2) as well as from lower bound in integrals over 
the terms m^/r^, or 1/s. In these latter cases the dependence on niq enters F2{x, Q"^) through 
the multiplicative factor 

_coll 

c = l- ™" = 1 S. (Q 4) 

M2 M2(l-x)' 

where we can set ruq = in (6.4). Moreover, as the ratio t^SiI/-^^ formally of power correction 
type, we can neglect it even for niq > 0. The term ln(l — x) in c!f^ causes problems in 
phenomenological parameterizations primarily because it appears there decoupled from the value 
of the quark mass nig with which it originally entered the expression (6.2) for Tmin and thus 
persists there even in the limit rUq 0, or when the quark mass is replaced by the initial scale 
Mo. 

In the presence of the lower cutt-off Mq on the integration over the quark virtualities r^JJ 
in (6.2) should perhaps be replaced with tq = max(my/(l — x), Mq). Sending rUq 0, but even 
for all realistic values of light quark masses and accessible x, we get tq = Mq. This implies the 
replacement ln((l — x)/x) —>■ ln(l/a;) in c!y^\ thereby removing most of the practical problems 
with the term ln(l — x) in (2.17). 

Beside ln(l — x) there is another term in Cy°^ that comes from nonzero quark mass, namely 
the last "+1" in the nonlogarithmic part 8x(l — x) — l = 8x{l — x) — 2 + 1. Discarding also this 
remaining trace of rUq we get, setting = Q^, 

C^°) (x, 1) = [x^ + (1 - xf] In ^ + 8x(l - x) - 2. (6.5) 

This expression is close to, but not identical QED formula for the constant part of F2{x,Q'^) 
obtained for massless quarks coupled photon with nonzero virtuality > [24] 

C^°) (x, 1) = [x^ + (1 - xf] In ^ + 8x(l - x) - 2. (6.6) 

The only difference between (6.5) and (6.6), i.e. the additional 1/x in the logarithm \n{\/x'^) 
with respect to (6.5) reflects the fact that for niq = and P^ > the lower bound r^jJJ = xP^. 
The expression (6.5) for C^^ is also close to that used for different reasons by Schuler and 
Sjostrand in their SaSlM and SaS2M parameterizations 

C^°) (x, 1) = [x^ + (1 - x)^] In ^ + 6x(l - x) - 1. (6.7) 

In Fig. 6 all three expressions (6.5-6.7) are compared to the standard one of eq. (2.17). 



7. Conclusions 

I have presented a reformulation of LO and NLO QCD analysis of which separates genuine 
QCD effects from those due to pure QED and satisfies the basic requirement of factorization 
scale invariance. This reformulation differs substantially from the conventional LO and NLO 
analysis of F2. 
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Compared to the standard formulation at the LO it requires the inclusion of four additional 
terms, proportional to c'^^ , k^^ , C^^ and C^"^ , all of which are known. Whereas in the standard 
approach the first three of them are part of the NLO approximation, the last one, i.e. the one 
proportional to enters the standard formulation first at the NNLO. Detailed numerical 

analysis of the contributions of all these terms shows that in most of the kinematical region 
accessible experimentally their sum is much more important than the contribution of the term 
included in the standard LO expression for F2 and induced by ■ It is shown that in the 
reformulated LO QCD analysis of F2 the part of C^*^^ proportional to ln(l — x), which in 
the standard formulation causes problems in the region x ^ 1, is much less troubling. All 
the quantitative considerations were carried out for the pointlike part of F2 in the nonsinglet 
channel and under the simplifying assumption /3i = 0. A realistic analysis of experimental data 
will necessitate the extension of the present formulation to the pointlike part of F2 in the singlet 
channel as well as the addition of the hadronic components in both channels. At the NLO a 
complete analysis requires the knowledge of two so far uncalculated quantities and is therefore 
currently impossible to perform. 
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